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FOREWORD 


This interim report was prepared for NASA Jet Propulsion Laboratory. 
Project Monitor on this contract is Richard M. Clayton, Liquid Propulsion 
Section of Jet Propulsion Laboratory, California Institute of Technology, 

Pa sadena , California . 

This effort was conducted for the Contract NAS 7-467 for the period 
May 1, 1966 to May 30, 1967. Dynamic Science number assigned to this 
report is SN-87 . 

This is an interim report (12 months) covering a continuing research 
program being conducted to develop a meaningful model of random wave 
phenomena occurring with hypergolic propellants during starting transients 
and during steady operation. 
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I. SUMMARY 


High amplitude pressure waves (pressure spikes) associated with 
vacuum engine starting transients, and pressure pulses (pops) occurring during 
steady state operation of hypergolic liquid rocket engines can be totally or 
partly detrimental to the rocket engine. The purpose of this research is to 
describe the conditions which, if present in a rocket combustion chamber during 
starting transient and during steady operation, will produce the observed pressure 
disturbances. Thus, the mechanism by means of which pressure spikes and pops 
develop are studied, however, the definition and verification of the pressure 
spike mechanism is emphasized. Both theoretical studies and laboratory 
experiments are utilized to accomplish this end. 

A detailed model is being developed which interrelates equations involv- 
ing: 1) transient flow rate, 2) chamber pressurization , 3) preignition chemistry 
and delay, and 4) detonation processes. The model will determine conditions 
which exist in the chamber until ignition and, ideally, will quantitatively predict 
the magnitude and time of occurrence of the overpressure. 

To date, certain of the individual system processes of the model have 
been formulated and their solution has been obtained with the aid of a digital 
computer, whereas other processes still require additional analysis and experi- 
mentation. 
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II. INTRODUCTION 


The importance of gaining a better understanding of starting transients 
and steady-state operation in hypergolic liquid rocket systems is borne out 
by the many applications of hypergolic propellant systems. 

A troublesome, and at the same time, highly complex phenomenon which 
is observed during engine starting transients is the occurrence of high pressure 
waves, usually termed pressure spikes. Similarly, pressure dlstrubances , 
commonly referred to as pops, are known to occur during steady-state rocket 
engine operation. Various attempts to explain the nature of these pressure 
disturbances have been made; however, most of these attempts are qualitative. 

The main purpose of this investigation is to develop an analytical model 
describing combustion chamber conditions during the start transient, although 
chamber conditions leading to pops are also studied. By necessity, the analyti- 
cal model is semi-empirical in nature, that is, containing several parameters 
which are best evaluated from experiment, or empirical correlation. 

In order to explain these pressure disturbances , it has generally been 
assumed that explosions of accumulated propellant masses are the cause of the 
observed pressure waves. As the result of more recent experimental studies, it 
was found that under certain conditions, detonable chemical reaction inter- 
mediates accumulate in the combustion chamber. This fact, evidently, adds a 
further complicating aspect to the quantitative description of pressure spikes. 

It would be clearly presumptuous to claim that a model based on theory 
alone could provide results of highly practical value to designers in order to 
eliminate or control spiking. Thus, the intent and contribution of the present 
study is to provide an analytical framework which will make possible parametric 
studies, showing which parameters affect chamber behavior most significantly 
and thus, will suggest the measurement of these parameters. 

The present model is one -dimensional and is based on time dependent 
differential equations formulating physical and chemical processes 
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governing chamber conditions. The sets of equations consist of: 

(1) Transient flow equations with time dependent resistance, 

in prf^nnp anH f 1 ow ppnaplfanpp 

(2) Vaporization, freezing, and pressurization equations. 

(3) Chemical kinetic equations including concentration and 
temperature dependence. 

(4) Detonation equations for a multiphase system. 

The actual solution of these equations is carried out numerically by 
a finite difference method and the computer program is arranged such that each 
of the cited sets of equations is treated as one unit. It is clear that the 
hydrodynamic and physiochemical equations are interdependent. For example, 
calculation of chamber pressure requires knowledge of the injection flow rate and 
the entering liquid and vapor mass fraction, temperature, preignition chemistry 
with related gas production, propellant vaporization rate, and nozzle exhaust 
losses. Changes in chamber pressure, in turn, influence the injection flow 
rate. 

In the ensuing discussion the problem of pressure spikes and pops will 
be defined and a background will be provided in terms of a literature review 
and some order of magnitude calculations . The ignition spiking analysis will 
be discussed in detail and the results of certain parametric computational 
studies will be presented. 
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III. PROBLEM DEFINITION 


Discussion 

The pressure spike phenomenon is characterized by a sudden rise in 
combustion chamber pressure during the ignition process in a rocket engine. 
Pressure rise rates of the order of 200 psia//i sec are observed, and peak 
pressures as high as 5000 psia have been recorded within the combustion 
chamber. Such behavior, represented schematically in the sketch below, 
is especially likely to occur at vacuum engine ignition. 



Although a considerable variance of the pressure spike amplitude is observed, 
the duration of a pressure spike is consistently of the order of one millisecond. 
The overall pressure spike has been resolved experimentally into a fine 
structure, namely, into several separate pressure peaks, each occurring 
over a duration of 70 to 100 /isec. 
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In contrast to the high amplitude (potentially destructive) pressure 
spikes associated with starting transients, the relatively low level pops 
appear to have been studied to a lesser degree. The characteristics of a 
chamber pressure pop resemble those connected with nonlinear, pulse or bomb 
induced disturbances. For example, the similarities consist of the steep- 
fronted pressure rise, the resulting maximum overpressure, and the damping 
characteristics. Both structural -hydrodynamic and physical-chemical effects 
have been assumed to have a part in popping. Experimentally it has been found 
that popping is influenced by the propellant combination, mixture ratio, and 
manufacturing peculiarities of a specific injector. 

In general, any propellant system can exhibit overpressures at particular 
operating conditions. In the following table are listed hypergolic systems which 
are known to exhibit transient overpressures. 


TABLE I. HYPERGOLIC PROPELLANT SYSTEMS EXHIBITING 
SPONTANEOUS PRESSURE DISTURBANCES 


Propellant Systems 

Spiking 

Popping 

1. Nitrogen tetroxide - Aerozine 50 

X 

X 

2. Nitrogen tetroxide - Monomethylhydrazine 

X 

X 

3. Nitrogen tetroxide - Hydrazine 

X 


4. Nitrogen tetroxide - Unsymmetrical 
Dimethylhydrazine 

X 


5. Nitric acid - Ethylenediamine 

X 


6. Nitric acid - Triethylamine 

X 


7. Nitric acid - Hydrazine 

X 


8. Nitric acid - Unsymmetrical 
D im e thy lhydra z in e 

X 


9. Nitric acid - Dimethylaminoethanol 

X 


10. Nitrogen tetroxide - 50% hydrazine - 

50% monomethylhydrazine 
1 

X 

4 

X 

J -1 
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These chemical systems are known to form unstable intermediate 
compounds under low temperature--low pressure conditions. It is possible 
that these intermediate compounds play an important role in the spiking 
phenomenon. 

The literature dealing with overpressure phenomena in the cited hyper- 
golic propellant systems is quite extensive. Rather than provide a complete 
listing of references treating the problem, only the most pertinent and most 
illustrative will be presented. The aim is to supplement the problem definition 
by listing and commenting on existing experimental work and existing hypotheses. 

Literature Review 

Reference 1. - R. M. Knox, S. J. Minton, and E. B. Zwick, Space 
Ignition, AIAA Second Propulsion Joint Specialist Conference, Colorado Springs, 
Colorado, June 13-17, 1966. 

Knox, Minton, and Zwick report very recent information on pressure 
spikes. They group pressure spikes under low backpressure conditions into 
three groups: 

(1) Vapor phase detonations with resulting peak pressures of less 
than 1500 psi. Frequently, smooth ignition occurs without an 
indication of detonation. No case of engine damage resulting 
from this type of ignition has ever been observed at Marquardt 
(no theoretical calculations are given to support the fact that 

a vapor phase detonation occurs.) 

(2) Large pressure which can damage an engine has been traced to 
a condensed phase explosion. This type of explosion is be- 
lieved to be due to the accumulation of frozen fuel on the wall 
of the chamber. Thus, if at the start of the next pulse the 
oxidizer valve is opened before the fuel valve, the oxidizer 

is believed to form explosive compounds with the frozen fuel 
in the engine. Residual frozen fuel can be present when the 
engine is started as a result of emptying from the dribble 
volumes, as a result of a leaky valve, and as a result of 
improper test procedures. 

(3) Propellant manifold explosions - which do not occur under 
space conditions. 
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Reference 2 . - S. J. Minton and E. B. Zwick, Hypergolic Combustion 


Initiated at Low Pressure, Aviation Space Conference of the American Society 
of Mechanical Engineers, March 1965, 47 pp. 

Reference 3 . - J. J. Kappl and R. M. Knox, Altitude Ignition of 
Hypergolic Bipropellant Rockets , Marquardt Corporation, 1965, 40 pp. 

These references present information obtained in the Marquardt 
Corporation's program to investigate the nature and prevention of high pressure 
spikes during ignition of 50% hydrazine-50% UDMH (Aerozine 50) and nitrogen 
tetroxide fueled rockets at low ambient pressures. Significant experimental 
results are: 

(1) The magnitudes of the peak pressure follow a random distribution. 

(2) The peak pressure produced seems to bear no relationship to the 
time delay but is randomly distributed. 

(3) The highest peak pressures occur with a 4 millisecond oxidizer 
lead. 

(4) The effect of increasing the propellant temperature is to shift 
the distribution of pressure peaks to lower values. 

Photographs of the ignition process are presented which are claimed to 
be suggestive of an initial deflagration followed by a detonation of the remaining 
material . 

Reference 4 . - M. L. J. Bernard and J. Dufour, On the Existence of 
Detonation Conditions in the Combustion of Some Nitric Acid Propellants, 

Eighth International Symposium on Combustion, 1960, pp. 1074-1084. 

Bernard and Dufour conclude that in the nitric acid-fuel (amines, 
alcohols, hydrazines, etc.) system there is a close relationship between the 
formation of unstable intermediate compounds and pressure spikes. Various 
explosive nitrates were found on the chamber walls of the test apparatus used. 

In the nitric acid-furfuyl alcohol system the rate of pressure rise for a normal 
ignition is 17.8 Kg/cm 2 /msec or 252 psi/msec. The temperature rise was 
estimated at 100°C over the entire ignition delay period. These pressure and 
temperature rises are believed (by Bernard and Dufour) sufficient to cause 
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deflagration and detonation respectively, of an explosive nitrate. A Fastax 
film study showed that an "abnormal" ignition (pressure spike) occurred when 
the chamber was first filled with vapors causing the entire chamber to light up 
at once. Contrasted to this was a "normal" ignition which occurred at several 
points in the chamber at the same time. (This suggests that perhaps a vapor 
phase detonation is initiated by a local hot spot or perhaps an explosive 
intermediate compound.) 

Reference 5 . - R. E. Martens, Investigation of the Hypergolic Ignition 
Spike Phenomena, McDonnell Aircraft Corporation, Internal Report, 27 pp. 

Martens carried out an experimental investigation of vacuum firings 
of a 25 lb. thrust engine using nitrogen tetroxide and monomethylhydrazine . 

The preignition accumulation of liquid oxidizer on the thrust chamber walls, 
verified experimentally, is claimed to be the source of the severe (^ 400 psia) 
ignition spike phenomena. Upon fuel injection an overspray collects on this 
coating in a small segment of the chamber. Following a vapor phase ignition, 
originating in the uncollected propellants, these accumulated propellants react 
providing a locally high pressure. This high pressure region, bounded on both 
sides by a lower pressure, generates a wave system analogous to that of a dual 
diaphragm closed shock tube. 

A simplified ignition spike analytical model incorporating the experimental 
findings was developed. The model was used to calculate the relationship be- 
tween maximum ignition spike pressures and oxidizer lead time which duplicated 
the experimental results. 

Severe spikes were only encountered during oxidizer lead starts at 
vacuum conditions and their severity increased with decreasing propellant and 
hardware temperature. 

In the McDonnell injector configuration the orientation of the oxidizer 
orifices was such that the oxidizer streams were canted outwards. Only an 
oxidizer to fuel ratio of 1 . 6 was used in the tests. 

On the fuel lead side severe spikes were not encountered, only gradual 
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pressure rises to approximately 400 psia were measured. On the oxidizer 
lead side, the ignition spike level was found to increase with lead time until 
a plateau was reached. This plateau, which persists for approximately 14 
to 15 milliseconds, is then terminated by a regime in which ignition spikes 
could not be generated. 

Reference 6 . - Norm Chaffee, NASA Manned Spacecraft Center, 

Houston, Texas, verbal communication. 

Mr. Chaffee reports that the Manned Spacecraft Center is carrying out 
an experimental investigation of the spiking problem using the Marquardt 
engine. He mentions that a condensed phase residue is present in the chamber 
after firing and that hydrazinium nitrate (N H NO ) was identified. 

^ D j 

Reference 7 . - Thiokol Chemical Corporation, Reaction Motors Division, 
Hypergolic Ignition at Reduced Pressures, Technical Report No. AFRPL-TR-64-175 , 
December 1964. 

Reaction Motors Division investigated the effect of several parameters 
on the ignition delay of propellant combinations containing ^^4 and IRFNA 
as the oxidizers and unsymmetrical dimethyl hydrazine, monomethyl hydrazine 
and 50% UDMH-50% h^H^ as the fuels. Parameters investigated were injection 
velocity, impingement length and angle, manifold feed configuration, propellant 
temperature, and ambient pressure. It was found that ambient pressure alone 
had the only significant effect on the ignition delay characteristics of all the 
injector and environmental parameters investigated in unconfined impingement 
tests . 

Reference 8 . - B. R. Lawver, and J. J. Kappl, Effect of Additives on 
Altitude Hypergolic Ignition, AIAA Second Propulsion Joint Specialist Conference, 
Colorado Springs, Colorado, June 13-17, 19 66, 31 pp. 

This work experimentally evaluated the effectiveness of selected fuel 
additives in eliminating or reducing altitude ignition pressure spikes in the 
propellant combination ^O^/Aerozine 50. Eleven diffexeuL aduitiv cS were mixed 
with the fuel in the proportion of 1% by weight. None of the chemical additives 
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tested were effective in reducing the maximum pressure spike experienced with 
^C^/Aerozine ^ * However, three of the additives significantly increased 
both the mean and the maximum pressure spike level, thus indicating that addi- 
tives can influence the altitude ignition process. 

Reference 9 . - B. R. Lawver, Effect of Ignition O/F on Spike Pressure 
Distribution, Marquardt Corporation Interoffice Memo 153-75/49, May 1966. 

Results are reported of spike tests conducted with cavitating venturi 
tubes placed in the propellant valves to shift the O/F at ignition. Certain 
off O/F ratios resulted in reduced ignition spike pressures. 

Reference 10 . - G. B. Skinner, W. H. Hedley, and A. D. Snyder, 
Mechanism and Chemical Inhibition of the Hydrazine-Nitrogen Tetroxide 
Reaction, ASD-TDR-62-1041 , December, 1962. 

Reference 11 . - H. G. Weiss, A Basic Study of the Nitrogen Tetroxide - 
Hydrazine Reaction, Dynamic Science, SN-4500, July 19 65. 

These reports are representative of several that mention the formation 
of intermediate compounds in the nitrogen tetroxide-hydrazine reaction. Weiss 
reports the formation of hydrazinium nitrate (N^^NO^) in a low temperature 
reaction between liquid nitrogen tetroxide and liquid hydrazine. Skinner reports 
the formation of ammonium nitrate (NH 4 NC> 3 ) in a preignition reaction between 

dilute vapors of N.O , and N_H . . 

2 4 2 4 

Reference 12 . - R. Friedman, W. P. Barnes, and M. Markels , Jr. 

A Study of Explosions Induced by Contact of Hydrazine - Type Fuels with 
Nitrogen Tetroxide, ASD-TDR-62-685 , September 1962 . 

Friedman, Barnes, and Markels studied the occasional explosions which 
result when a hydrazine -type liquid is brought into contact with liquid nitrogen 
tetroxide with little confinement. The authors propose the existence of an 
adduct that forms on propellant contact and would be responsible for the 
explosive reaction. This type of explosion may be related to the pressure spikes 
observed with this propellant system. 
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Reference 13 . - Henry E. Perlee and Theodore Christos, Summary of 
Hypergolic Ignition Spike Phenomena, U. S. Department of the Interior, 

Bureau of Mines Final Report No. 3982, April 8, to December 31, 1965. 

The Bureau of Mines carried out a literature survey of hypergolic 
ignition spike phenomena for the MSC in Houston, Texas. One hundred and 
sixty three articles are reviewed including some of these mentioned above. 

Topics covered include: physics of spray formation, hydrodynamics of fluid 
streams, jets in initially low pressure environments, volatile liquid flow in 
initially low pressure piping, heterogeneous pipe flow, hydrodynamic instability 
criteria of liquid jets, single particle physics, cooperative phenomena, drop 
kinematics, transport phenomena, physics of sprays, chemistry of combustion, 
ignition, flame propagation, and gas dynamics. It would serve no useful 
purpose to comment on all of these articles here. Most of these refer to physical 
and chemical processes which take place in an engine and which could have 
relevance to the hard start phenomena. 

Literature on popping is rather scarce. The following reports give 
experimental results on popping: 

Reference 22 . - Richard R. Weiss and Raymond D. Klopotek, 

Experimental Evaluation of the Titan III Transtage Engine Combustion Stability 
Characteristics, AFRPL-TR-66-51 , March, 19 66. 

Reference 23 . - R. R. Weiss, T.J.C. Chew, and Lt. R. D. Klopotek, 

A Combustion Stability Evaluation of Various Hydrazine and Hydrazine Blend 
Fuels, Air Force Rocket Propulsion Laboratory, Edwards, California, 42 pp. 

This work examined the propellant combinations N o 0.-A50, N„0,-MMH 

2 4 2 4 

and N^O^-Monozine 50. The significant results of the work which involve 
popping are: 

(1) There appears to be no correlation between baffle cracks and 
the occurrence of popping. 

(2) Experimental results did not support the speculation that large 
amounts of film coolant are a major cause of pops. 
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(3) N^O -A50 and appear equally susceptible to the 

popping type of instability. The N^O^/Monozine 50 combina- 
tion showed a considerable increase in susceptibility to 
popping. 

Reference 24 . - R. M. Clayton, Combustion Roughness and Dynamic 
Stability, Jet Propulsion Laboratory, a section of a future report. 

Clayton shows that popping is related to conditions in proximity of 
the chamber walls . The absence of reactants throughout a substantial space 
near the chamber boundary and especially in the early reaction zone near the 
injector was the primary contributor to the absence of popping in a particular 
engine configuration. 
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Order of Magnitude Calculations 

In order to evaluate the overall relative importance of one-and two- 
phase detonation and explosion, order of magnitude calculations were carried 
out applicable to spiking and popping. In order to make these calculations tract- 
able certain simplifying assumptions had to be made. For example, only one- 
dimensional aspects of detonation and explosion are considered, and complete 
reaction of the propellants is assumed in the case of two~phase detonation. 

Order of Magnitude Calculations Applicable to Spiking . - 

Several order of magnitude calculations were made using Marquardt 
rocket engine data (Refs. 2, 3, 8, and 9). 

Gaseous Detonation. The pressure ratio, P /P , was determined 
for the case of a vapor phase detonation. The NASA Lewis Chapman-Jouget 
detonation program was used. This program utilizes a Newton-Raphson 
iteration procedure to calculate the pressure ratio across the detonation wave. 
Initial conditions were taken as: initial temperature, 278°K, initial pressures, 

2 psia and 4 psia. The resulting pressure ratios were: for 2 psia initial 
pressure, 7 2.2; for 4 psia initial pressure, 73.8, based on O/F = 2. 

The significance of these results is that these pressure ratios are 
not high enough to account for pressures obtained in the Marquardt engine 
(Ref. 9). 

Heterogeneous Detonation. Heterogeneous detonation calculations 
were made using the theory of Ragland, Dabora , and Nicholls (Ref. 14). 

Assuming that the propellant drops bum completely and that the Chapman-Jouget 
condition is realized, conservation of mass, momentum, and energy equations are 
obtained which lead to the following expressions: 




1 + oo'jMj 

1 + To 
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where subscript 1 indicates initial conditions, and 2 final conditions, and where 
a. = (p^+ P] )/ p 

M = Mach number 

p' = mass of liquid per unit volume of mixture 

Xj 

Q = heat release due to chemical reaction minus phase change 
y = ratio of specific heats 
e^ = internal energy of the liquid. 

Order of magnitude calculations were made giving as a function of 
O/F ratio. 

The NASA-Lewls Chemical Equilibrium program was used to determine 
the heat of reaction and equilibrium product gases y . The initial temperature, 

T^ , was taken as 40°F, the initial pressure before spiking as 2 psia (see Figure 
1) and y = 1.4. For this order of magnitude calculation the internal energy 
of the liquid and the heat change due to phase transition were neglected. The 
initial gas density was determined from the ideal gas law at the initial condi- 
tions . A quality of 0.043 lb vapor per lb propellant was used (Ref. 2). was 
determined as a function of O/F. Figure 1 compares the theoretical results 
with experimental results obtained at Marquardt, (Ref. 9). Since the pressure 
ratio, P^/P ^ , varies with quality. Figure 2 was prepared to show the final 
pressure as a function of quality with initial pressures of 2 psia, 4 psia, and 
8 psia as parameters, and an O/F = 2. 

In addition to the equations of Ragland, Dabora, and Nicholls (Ref. 14) 
the following expressions were used in obtaining the curve: 

M 

f = 

+M 
L v 
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M. + M 

l V 

P M 


V +v 

l V 

P M 

= 

h 

»-£)+£ p/ p v 

where: 



f 

— 

quality of vapor (pounds of vapor per pound of 
propellant) 

a 

= 

P v d- f )+ 

p v (1 -f ) +£p i 

M 

V 

- 

mass of vapor 

m l 

= 

mass of liquid 

m t 

= 

total mass of propellants 

V T 

= 

total volume of propellants 

V i 

= 

volume of liquid 

V 

V 

= 

volume of vapor 


= 

density of liquid 

p v 

= 

density of vapor 

P M 

= 

density of mixture 

The results of these order of magnitude calculations, as shown in 


Figures 1 and 2 indicate that heterogeneous detonation can account for 
pressure spike amplitudes observed in the Marquardt engine. 

Accumulation of ProDellant. - A constant volume combustion calculation 
can be made for the Marquardt engine. Thus, PV = (W/MW)RT, where 


. W 
R 
T 
V 

MW 

yields P 


_0 

0.965x10 (assuming a 2 millisecond residence time) 
1545 ft lbf/lbmole °F 

5850°R (from NASA-Lewis program with Q/F= 2) 

6.28 in. 3 

22 (approximate product gas molecular weight from 
NASA-Lewis program for 0/F=2) 

760 psia. 
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An absolute maximum pressure attainable by constant volume combustion may 
be estimated by assuming no nozzle exit flow throughout the ignition delay 
period. Thus, on the basis of the above conditions, 

A P/At = 380 psia/millisecond of ignition delay . 

For a 5 millisecond ignition delay time one obtains 1900 psia. 

For the case of no nozzle exit flow, a pressure is obtained which is 
in reasonable agreement with some observed Marquardt pressure spikes. How- 
ever, since the assumption of no nozzle exit flow is unrealistic, the calculated 
value of 7 60 psia for the final pressure attained by burning all of the propellant 
within a 2 millisecond residence time tends to rule out explosion as a mechanism 
for the hard start observed in the Marquardt engine. 

Quality of Marquardt Engine. A value for the quality (pound vapor per 
pound propellant) was determined from the following data: 

Residence time, based on 

injection velocity 1.5 to 2.0 millisecond 

Volume of chamber 6.28 in. 3 

_3 

Cold flow rate 0.483x10 lbm per millisecond 

3 

Thus, p^ (density of the mixture) varies between 0.212 lb/ft and 
0.266 lb/ft 3 . Using the relationship = p /£(l-f)+f o^/ cy],the quality (f) 
varies between 0.097 and 0.12 or an average quality of 10.8 per cent. 

Explosion. Using the ideal gas law relationship, P = p^RT, and the 
variation of the mixture density with quality, p^ = p^/[(l -f)+f p^/ , Figure 3 
was prepared to show the final pressure, P, as a function of quality, f. An 
explosion is here defined as instantaneous reaction yielding a pressure calcu- 
lable from the ideal gas law. A final temperature of 5850°R and an O/F of 2 
was used. Figure 3 shows that a relatively high liquid mass fraction is re- 
quired to yield explosion pressures of the same order of magnitude as are 
found experimentally. 
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Order of Magnitude Calculations Applicable to Popping. - 
Heterogeneous Detonation. Employing Dynamic Science's steady-state 
combustion program in conjunction with Ragland, Da bora , and Nicholl's detona- 
tion model (Ref. 14), the detonation pressure ratio was calculated as a function 
of distance from the injector face. For an O/F = 2 and the N^O^/MMH system, 
the fraction of propellant vaporized as a function of distance from the injector 
face was calculated using Dynamic Science's combustion program. These data 


are: 


and 


Distance from Face 
Inches 


Per cent 

Unbumt Material 


0.25 

0.50 

0.75 

1.00 

3.00 

5.00 

7.00 

9.00 

The equations: 






1+ay^M 
>■ + y n 


93.0 

88.0 
82.0 

78.0 
57.5 

47.0 

40.0 

36.0 


C _T. 
Pi 1 


were then used. Q was modified by taking the heat of reaction from the 
NASA-Lewis program and multiplying by the per cent of unbumed material. The 
value for was taken as 1.4. The specific heat ratio, y ^ , and the steady-state 
temperature, T^ , were taken from the NASA-Lewis Program. Figure 4 shows 
the large effect which percent unbumt material and droplet number concen- 
tration can have on detonation overpressure near the injector face. 

Explosion Phenomena. Besides the magnitude of peak pressures 
developed in an explosion, the characteristics of pressure decay are of significance 
On the basis of published pressure data available at Jet Propulsion Laboratory, 

various distances from a bomb exploded in the center of a rocket engine, it 
20 





is possible to compare the experimental with isentropic and isothermal expansion. 
Employing the ideal gas law and the following initial conditions: 
weight = 13.5 grains of RDX 
grain size = 0.279 in O. D. , 1 . 5 in. long 
75 

Tw 6000° 

one obtains a pressure of 114,000 psia at a radius of 0.49 inches (based on a 

sphere volume equivalent to the cylinder of explosive). Using this as the 

initial point,the experimental data are plotted in Figure 5 . Also plotted are 

data for an isothermal expansion (PV = constant) and isentropic expansion 

y 

(PV = constant). According to Figure 5, at radii near 4 inch or greater, the 
experimental pressure decay curve approximates isothermal expansion. 

The results presented in Figure 4 demonstrate the sensitivity of the heterogeneous 
detonation pressure ratio to distance from the injector. Pressure ratio resulting 
from a vapor phase detonation, on the other hand, does not vary significantly 
with distance from the injector. 

The significance of the order of magnitude calculations applicable to 
spiking and popping can be summed up in the following: 

(1) Gaseous detonation by itself does not account for peak pressure 
characteristic of spiking, 

(2) Heterogeneous detonation calculations indicate that spike pressures 
are attainable in regions where the liquid propellant mass fraction 
is large, 

(3) Although constant volume explosion for high liquid to vapor mass 
ratio yields pressures of spiking magnitude, the closed system 
assumption for the combustion chamber is unrealistic, and 
therefore, it is unlikely that explosion accounts for spiking, 

(4) Vapor phase and heterogeneous detonation pressures approach 
nearly the same value at relatively large distances from the 
injector face, that is, at locations where the liquid mass fraction 
is low. The magnitudes of these detonation pressures are con- 
sistent with overpressures characteristic of pops. 

(5) In general, the order of magnitude calculations carried out point to 
the need for a more detailed modeling of the overpressure phenomena. 


22 







The analytical framework to be described in the following section makes 
possible a more realistic study of the relative importance of various physical 
and chemical processes involved in hypergolic propellant ignition and 
combustion . 
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IV. IGNITION SPIKING ANALYSIS 


The pressure spiking phenomenon, as defined in this study, is closely 
associated with the engine ignition transient. Thus, it is important to consider 
the ignition sequence variables in order to understand better the complex inter- 
action of the physical and chemical processes. 

Discussion of Typical Starting Sequence 

Since vacuum or high altitude ignition creates, in practice, the most 

susceptible condition for spiking, such an ignition sequence (Ref. 8) will be 

employed to illustrate certain measurable processes taking place in a 100 lb. 

thrust altitude control engine. In referring to Figure 6, describing a typical 

altitude ignition sequence, attention should be called to the fact that the 

various valve signals may be either electrical or real injection lead signals. 

The actual injection lead, of course, is of particular interest in connection 

with the ignition delay. In Figure 6, the initial chamber pressure is the cell 

pressure (*w0.5 psia). According to the valve sequencing in Figure 6, the 

oxidizer leads the fuel by msec and about 7.5 msec after the oxidizer 

valve signal, the oxidizer starts to enter the chamber, and the pressure rises 

to about 1 psia. The fuel starts to enter the chamber about 7.5 msec after the 

fuel valve signal. In order to describe mathematically the altitude ignition 

process, it is first necessary to formulate a general framework which takes 

into account all important variables and their interrelationship. Figures 7, 8, 

9, 10, and 11 show schematically the complex interaction of the physical and 

chemical processes. The chamber pressure P (t) is the primary variable which 

o 

is being described before and during ignition. The chamber pressure is affected 
by the tank pressure through the four simultaneous and/or parallel systems 
shown in Figure 7. The four systems being described mathematically are: 

(1) The transient flow process (Figure 7), 

(2) The chamber pressurization process (Figure 6) 
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Fuel Valve 
Signal 



Time = msec 


FIGURE 6. Rocket Engine Pre-Ignition Pressure History Data 
Taken From Reference 8. 
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FIGURE 7. TRANSIENT FLOW SYSTEMS (I) 


PREIGNITION LOSS 



FIGURE 8. CHAMBER PRESSURIZATION SYSTEM (II) 
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FIGURE 9. PREIGNITION CHEMISTRY SYSTEM (III) 



FIGURE 10 


ROCKET CHAMBER IGNITION 
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FIGURE 11. COMBUSTION/DETONATION SYSTEM (IV) 



(3) The preignition chemistry and ignition delay (Figure 9), and 

(4) The ignition process (Figure 10). 

These four systems are described in this report along with a section 
on the chamber pressurization-ignition delay interrelationship. 

In the transient flow system of Figure 7, the tank pressure is the 
input while the vapor, V^(t), and liquid, L(t) flow rates into the chamber 
are the computed outputs. 

The chamber pressurization system receives an input liquid flow at L(t) 
and through the vaporization process shown in Figure 8, converts this 
liquid into vapor V 2 (t) and condensed phase material C(t). The complete 
pressurization system must also account for initial vapor pressure flow 
through the injector orifice and preignition chemistry gas generation. The output 
of this system is the preignition chamber pressure P c (t), and temperature T c (t). 

The preignition chemistry system of Figure 9 must be treated in 
greater detail than is available in the literature. Intermediate formation 
and gas generation contributes significantly to chamber pressurization. For 
example, with an oxidizer lead, high speed motion pictures indicate the 
presence of colorless 1^0^ gas prior to fuel injection. However, almost 
instantaneously with fuel Injection into the oxidizer vapor, redish brown NC^ 
appears in the vapor space. This phenomenon may be due to heat generation 
or intermediate gas formation, either of which affect the chamber pressurization 
system of Figure 8. Thus the output of the preignition system will be the 
specification of intermediate vapor formation, Vg(t), and intermediate pro- 
ducts , I (t) , at P (t) and T (t) . 

c c 

The combustion/detonation system of Figure 11 involves an ignition 

or detonation decision based upon the incoming conditions of mixture, 

chamber pressure and temperature. This decision must be based upon 

propellant accumulation, heat loss, kinetics, and detonation transition 

length. Once the decision is made, the final spike produced, P s (t), or 

the combustion overshoot, P (t), may be calculated by methods presented 

c 
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in the order of magnitude calculations section. 

The chamber ignition process is driven by a tank pressure signal, 

P , which passes through numerous complicated processes of varying 

importance and produces a detonation spike P (t) or an ignition chamber 

s 

pressure, P (t). 

c 


Transient Flow. - The transient flow rate into the chamber must be 
known as a function of time in order to describe chamber pressurization. 

The technique being applied involves analysis of experimental 
pressure-tap data and empirical specification of a lumped parameter system 
which is valid over a limited range of operating conditions. 

Feed Line Dynamics. The range of conditions being treated at present 
are those encountered with pulse motor start-up at space-vacuum. The 
sketch below shows the diagrammatic setup. 



The tank pressure, P^ (t) , which may be assumed constant, is trans- 
mitted through the line impedance Z^ and results in the pressure P^(t), 
measured at the valve inlet and shown on the oscillograph trace below 

for the oxidizer P_ , _ (t) and fuel P.,_(t). 

MO MF 
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TYPICAL OSCILLOGRAPH TRACE 



The valve opening time in this trace is approximately one millisecond 
while the inlet pressure P^q continues to drop for approximately eight 
milliseconds. 

The header pressure will equal the liquid vapor pressure previous 
to header filling. The time of header filling is determined from the volume 
of the header and the respective mass flow rates through the valve and 
injector. Once the header is filled, its pressure will rise to an equilibrium 
value. 

By comparing both the header pressure and the chamber pressure 
with the vapor pressure, the appropriate mass flow equations can be applied. 
Thus, for the valve: 


w = 

K 

> 

s 

1 

< 

f° rp H 

= P V 

w = 

K 

/ P - p 
v M H 

f° rP H 

* p \ 


where: 

Py = vapor pressure 

P = header pressure 

ri 
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K 

For the injector: 
W = 


a constant 


V* 




M 


ZM ; '_g_ f 1 *-' p - p !c, -x_ J /y - 1 

R ^y+lV H V P__ V.y+1 / 

xi 


w = K ' /p H - p v 


w = k/p h - p c 


p > p , p £ p 
H V* C V 


P H > P V' P C > P V 


The pressure upstream of the manifold valve will be determined by a 
lumped parameter system while the flow rate into the chamber is determined 
by one of five sets of manifold/header/chamber pressure criteria. TRW recently 
had success in such a description by describing the feed line dynamics, on 
valve opening, with the following lumped parameter analog (Ref. 15). 



w,. 


The equations describing flow in this system are derived for a mass balance: 
dP 

_ a 

C — = w - w 
a dt a p 

dP h 

C, — = w - w 
b dt p b 

and a momentum balance: 

dw 


= gA (P„ - PJ - R (w^) 2 
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where: 


w = section inlet flow rate, lb/sec 

3 

w = section instantaneous flow rate, lb/sec 
P 

w, = section outlet flow rate, Ib/sec 
b 

= section inlet pressure , psia 

P. = section outlet pressure, psia 
b 

C = feed line capacitance, defined as 

r _ £AL 
144/3 

C = effective feed line capacitance at station (a) 
a 

C. = effective feed line capacitance at station (b) 
b 

/? = propellant bulk modulus , lb/in. 2 

1//3 = propellant compressibility, in.^lb, where 

l7 Wr 3 - 20xl0 ‘ 6ln - Vlb 

0 4 = 5.31 x 10~ 6 in. 2 /lb 

p = density of oxidizer = 90.31 lb/ft 3 

ox 

p^, = density of fuel = 54.73 lb/ft 3 

A = cross-sectional feed line flow area, in. 2 

L = length of line per section, ft 

g = acceleration due to gravity = 32.2 ft/sec 2 

R = lumped feed line frictional loss coefficient, ft/lb 


These three equations involve the input tank pressure P , and the 

3 

four unknowns: Output manifold pressure P, , and the three flow rates w , 

d a 

w , w, . The output flow rate w, is related to the manifold pressure and 
p b b 

the downstream header pressure by 


w. 


P, - P„ = R 
b H v 


A 

v 
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where: A^(t) is the valve open area and is the driving force which actuates 
the system. 

R is a variable nonlinear valve resistance for single and two-phase 
v 

flow. The downstream pressure in the header determines the type of flow and 
resistance existing in the valve and injector orifice. It is found by accounting 
for the mass accumulation in the header by applying a series of flow regime 
criteria . 

A slightly more complex computer program is being set up at the Manned 
Spacecraft Center in Houston (Ref. 16) which gives the flow rate into the 
chamber both with and without a preigniter. The model lumps feed line, injector, 
and combustor modules. The model was developed by both General Electric 
(Ref. 17) and the Manned Spacecraft Center (Ref. 16). Programming is nearing 
completion and this program could be used if desirable. However , this program 
does not treat two-phase and cavitating orifice flow, as is being attempted 
at Dynamic Science. 

Cavitation. In order to determine the existence of two-phase flow 
the photographs of Figure 12 were recorded with a 10 microsecond flash unit 
for the N O /NO flow system at the indicated time intervals after valve 
opening. The lighted portion of the first photograph shows a simulated two- 
dimensional manifold sandwiched between two pieces of plexiglass. In the 
next five photographs the N 2 0 4 /NC> 2 vapor is the dark substance while the 
light penetrates through the liquid. Two-phase flow and cavitation are present 
during the entire filling sequence and even in the last photograph the injector 
orifice is fully cavitated. The transient flow equations are controlled during 
cavitating periods by the two phase flow resistance R^. This proportionality 
constant will have to be determined experimentally as a function of back pres- 
sure and a digital program can then be written to solve the transient flow problem. 

Chamber pressurization. - The following differential equations apply to 
the low pressure transient flow process when the liquid is at the same temp- 
erature as the vapor. 
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FIGURE 12. N-O . MANIFOLD FILLING SEQUENCE 




where: 


V 27rr^n0 


2 V l nor 


; j 


per cent mass in drops smaller than given radius r 
drop radius 

the geometric standard deviation 
the mass median drop radius 


Evaporation: 


J = " P - P " a A 

. VTd gTg J d . 2irRT . 


G i = pounds mass of drop size i evaporated 

t = time 

^VTd = va P° r pressure evaluated at the temperature of the drop 
P = gas pressure 


Condensation: 


accommodation coefficient 

area of the drops 

molecular weight 

gravitational constant 

universal gas constant 

temperature of the drop 
dG, 


" = [ 


P - P 
gTg VTw 


c V 2n RT 


where: 


= mass of drops of size i condensing 

PyTw= vapor pressure of the substance at the wall temperature 

A = chamber wall area 

c 

T g = gas temperature 
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Mass efflux out nozzle: 


where: 

m 

7 

A* 



mass flow out the nozzle 


ratio of heat capacities, C /C 

P v 


the throat area 


Heat balance on liquid drops: 


TdG i = p e!’ 7r f 


C .dT.+ — irp.r 3 ,C dT, 
p-6 i 3 l i N-l ps is 


where: 

r = heat of vaporization 

G = pounds mass of drop size i evaporated 

p. = liquid density 

** th 

r = radius of the i drop group 

C ^ = heat capacity of liquid 

T. = temperature 

C ps = heat capacity of solid particles 
T\ = temperature of 1 th size solid drop 

Since P , P , and T, vary with time, it is necessary to use a 
v g d 

numerical integration procedure. Agosta (Ref. 18) uses such a procedure: 

Drop groups with specified mass distribution are determined from the 
logarithmiconormal distribution of Priem and Heidmann (Ref. 19). 

Drop groups (i.e. n drop groups) with a given mass distribution are 
chosen. The number of drops N. in each group is determined by 

N. = W ./ (r/3) rrr 3 p x 


whereiW^ = the mass of drops in each group. 

Based on Knudson-Langmiur kinetic theory, it can be shown that 
the rate at which ideal gas molecules strike a surface is given by 
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G - P 7 Mg/2trRT 

where : 

G = mass of molecules striking a surface per unit area per 
unit time 

P = gas pressure . 

R = gas constant 

M * molecular weight 

T * gas temperature 

g = gravitational constant 

If the surface is a liquid or solid condensate in equilibrium with 
its saturated vapor, and if the presence of the saturated vapor is denoted 
by P v then the above equation becomes 

G - P v 7 Mg/2ffRT 

where: P v * vapor pressure 

If a fraction a of the molecules condenses on the surface, the 
equation becomes 

G = a P v 7 Mg/2 ir RT 

where: a =» the accommodation coefficient 

The equation gives the rate of evaporation from a surface in contact 
with its saturated vapor when both surface and vapor are at temperature T. 

If some of the molecules that come off the condensate are not 
immediately removed so that a vapor cloud develops above the condensate 
surface, the net rate of evaporation is given by 

G = a(P -P ) 7Mg/2irRT 

v a 

where: P = the pressure of the vapor cloud from a single drop. Thus 
the mass of liquid evaporating in a finite time interval is: 



or 


G. . = 
ij 


G.. = 
ij 


[Wij N-l _P g^ T g^N-l J a ^ij^N-l A Vr(t !) 9 

d ij N-l 

r , W. At 

C [ P (T ) - P (T ) I 1 J N 

i|_ VVijN-1 g u g ; N-l 


where: = 

P v (T ) 

V ' i d ijN-1 

P (T ) 
g g N-l 


ij 

4t N = 

w. = 

J 


9 ^ ijN-1 


th 


a constant depending on the number of drops Initially in an i 
the vapor pressure measured at the temperature of the i tft size 
drop of the j group at the end of the N-l time interval. 

gas pressure at the end of the N-l time interval 

th th 

radius of the i size drop in the j group 


time interval chosen 


.th 


mass of liquid propellant injected during the j interval of 
time 

(T ) th th 

' d ijN-1 = the temperature of the i size drop of the J group at the 

end of the N-l time interval 

j = drop groups, 1 s] sn . j can start at a number greater than 
1 if the initial drops become completely evaporated 

th 

The total evaporation occurring during the N time interval is: 

n N 

(G ) = S E g 

v EVAPN ; N 3=1 pi ij 

At the end of every time interval, new drop radii and drop temperatures 
are calculated from: 


G ij 3 7Tp Ji L (r ij 3) N-l " (r ij ^N j 

and from the previous expression for G , , 

ij — , . ■*•/ O 

(r ) - (r , ! * JWllN-1 -Wn- 1 ^ At N ' 

1)N 1)N -1 p. r Vi 

ijN-1 diJN-1 


where: K 


a constant. 


group 
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The new droplet temperatures are determined from a heat balance: 


G ij T = ^ 3 ) np * f3 


(r 


and 


A ijN-1 pA dijN-1 
Kt [p_ /m t - P 


- t ) 
dljN' 


T dijN T dijN-l 


^w-r r ^N-J 


,r 


p i C pi r ijN-l J T dijN-l 


If liquid freezing occurs, the term x (4/3) it p ur 3 .,.,, . r fusion 

i ijN-1 

(where: x = the fraction of the liquid that has solidified) 

must be added to the right hand side of the heat balance equation. 

If the liquid completely freezes and falls below the freezing 
temperature, the term: 

(4/3) it P a r i j N _ 1 c ps ( T sl j N _ 1 - T si j N ) 

must also be added to the right hand side of the heat balance equation. 

A separate (G„, ._..)„ can be written for both the oxidizer and the 
EVAPN N ^ 

fuel. The mass of vapor condensed on the chamber wall during the N 
time interval would be: 


/ Mg 


(G COND ) N Pv (T s ) ] ttA c At N V2irR(T g ) N _ 1 

Nozzle Flow: For isentropic flow of a perfect gas through a nozzle with 

constant values of C , M, and y , the mass through the nozzle during the 

th p 

N time interval is: 


(M N° z) N = = 

u g'N-l 


***- &)' fl/7 " 1 


R 


th 


A mass balance on the vapor gives at the end of the N time 


interval 


(G EVAPN^N 


(G COND^N " * M NOZ*N (M c^N ~ * M c*N-l 
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where: 


(M ) = mass of vapor in the chamber at the end of the N** 1 time 

interval as determined from the ideal gas law 
V MP 

with V being the chamber volume. 

r* 


R T 


gN 


(ML t = mass of vaDor in the chamber at the end of the N-l time 
cN-l V C M P aN _! 

interval = — — — 3 

R T . . . 
gN-1 

In order to determine from the mass balance equation for the vapor, 

Agosta (Ref. 18) determines T XT by means of a mass weighted average: 

N n 


J=1 1=1 


G..(T.) .... . 
ij d ijN-1 


gN 
+ (MJ 


c'N-l T gN-l ~ ^ G COND^N T gN-l ^ M NOZ* T gN-l 


gN 


^ =1 G ij + ^Vn-1 " ^CcWn ^ M NOZ^N 


Ignition delay and preignition chemistry. - In order to determine 
the conditions which exist in a rocket engine at and prior to ignition, it is 
necessary to know the ignition delay time and the preignition chemistry. 

Knowing the pressure rise as a function of time, the ignition delay 
can be determined from a knowledge of the variation of ignition delay with 
ambient pressure. This information can be determined experimentally and in 
some cases is known. 

Reaction Motors Division (Ref. 7) has experimental data for the 
system ^O^-UDMH which fits the equation 

4n t d = 15.5 - 2.6 Zn (Pa) 

where 

T" d = ignition delay 

Pa = ambient pressure 

Although this particular correlation applies only to N 2 C> 4 -UDMH ( experimental 
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data of ignition delay versus ambient pressure is available for N o O,-A50, 

Z 4 

N 2 0 4 ~MMH and IRFNA-UDMH. 

A theoretical determination of the ignition delay can be based on 
the criterion that ignition occurs when the heat generated just exceeds the 
heat loss. This criterion is true for thermal ignition mechanisms, 
i.e. , no chain branching, and is probably true for hypergolic propellant 
systems . 

The rate of heat generation in a rocket engine can be expressed by: 


where 

AH 

V 




= AHVr + h .A (T -T) 
wl w w 

= the heat of reaction 
= volume of thrust chamber 

= Arrhenius reaction rate 

= heat of transfer coefficient from the chamber 
wall to the propellants 

= chamber wall area 

= wall temperature 

= propellant temperature 


The reaction rate is determined from: 


K C n C™ d -E//RT 
o F 


where: 


K 

C 


n and m 


E 


a constant 

oxidizer concentration 
fuel concentration 

experimentally determined exponents 
activation energy 


The heat loss can 


be determined from: 


= (hA + WC ) (T-T )+h A (T-T ) 
p' O' WZ w ' w' 
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a. 

‘loss 


where : 


A 

W 

C 


w2 


w 


w 


= heat transfer coefficient between gas in the thrust 
chamber and liquid and solid particles 

= combined area of the liquid and solid particles 

= flow rate out the nozzle 

= specific heat of reactant mixture 

= propellant temperature 

= temperature of incoming propellant 

= heat transfer coefficient from propellant to wall 

= chamber wall area 

= wall temperature 


The rate of change of temperature with time, based on an energy balance 
is. 


dT 

dt 


AHKC n C%- E/RT h .A (T -T) (hA+WC )(T-T ) 
o F wl w w p o 

" T " “ 


P C 


P 


m C 


m C 


h A (T-T ) 
w2 w w 

m C 


where: m = the mass of the system 

p = density of the system 

That portion of the heat generation term that deals with preignition 
chemistry (in this case was examined analytically with a pre- 

liminary kinetic model. 

The reaction rate of the low pressure, N^O^-N^H^ gas reaction 
was determined from the collision theory which states that: 

r ’ K °A C B 

-E/RT 


K 

A 

where K 


M ^ m a+ m b 
AB v 8?rRT M. M 


1/2 


= CT 2 


A *"*B 

= reaction rate constant, cc/ (molecule) (sec.) 


= concentration of nitrogen tetroxide, 
molecules/cc 
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A 

E 


- C „ - concentration of hydrazine, molecules/cc 
2 4 

= frequency factor 
= activation energy 


a 


AB 


= effective diameter of A and B upon collision 


CT „ +CT 

A B 


M , M = molecular weights 

A D 

R = gas constant = KgN^the product of Boltzmann's 

constant and Avogadro's number, ergs/(°K)(g mole) 

In this analysis a one to one molar ratio of reactants is assumed along 
with a reaction which is first order with respect to each reactant. 

For use in the equation, the collision diameters were calculated 
from the cube root of the respective molecular volumes. The molecular 
volumes of nitrogen tetroxide and hydrazine were determined from the 
sum of the respective atomic volumes (Ref. 20). Thus, 


CT N 2 H 4 = 4.25 x 10“ 8 cm = 4.25A 

a N o 0. = 4.7 x 10“ 8 cm = 4.7A 

2 4 

The activation energy, E, for the reaction was determined from 
one fourth the sum of the bond energies broken. In this case the sum of the 
bond energies broken is determined by postulating a reaction initiated by 
the following decomposition reactions: 

N 2° 4 (3) -*N 2 + 20 2 AH 2gea5 o K =-2.35Kcal/mole 

3N 2 H 4 (g) - 4NH 3 + N 2 -37 Kcal/mole 

Liquid hydrazine thermally decomposes to NH^ and N 2 at high 
temperatures in very short time intervals, such as transient times in 
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rocket engines (Ref. 21). The same decomposition reaction is assumed to 
hold for gaseous hydrazine. The activation energy is taken as one fourth 
the sum of the respective heats of dissociation of the above reactions, thus 
E = 9.84 Kcal/mole. 

Substituting in the above expression for reaction rate, one obtains: 

-a -8 1/2 --SS42 

-4.25x10 +4.7x10 1 , M , n n w in 7* 92+32 "] 1.98T- 

r l; 2 j ^8 (3. 14) (8. 3) (10 ) J e 

_ Molecules 
x C-„ TT } r — , or 

N 2 H 4 cc (sec) 

1 1/2 _ .10 -4969. 7/T _ _ moles 

r=1 ‘ 128T 10 e C N 2 0 4 C N 2 H 4 li^(^) 

If complete equilibrium reaction of takes place, water 

and nitrogen are formed as products: 


N 2°4 +2N 2 H 4" 4H 2 0 + 3N 2* 

The heat of reaction for liquid nitrogen tetroxide with liquid hydrazine 
is -248 Kcal/mole. However, experiments at Dynamic Science have shown 
that for liquid nitrogen tetroxide and liquid hydrazine reacting at room tempera- 
ture, the heat of reaction is of the order of 30Kcal/mole. Assuming this 
same heat of reaction for a gaseous reaction at room temperature, the 


rate of temperature change becomes: 


(30)I.128T^10 l0 C,, C^, . 


-49 69 . 7/T 
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Chamber Pressurization - Ignition Delay Interrelationship: 
The mass weighted average temperature 


^ G ij ^ T d^ijN-l + ^ M c^N-l^ T g^N-l ^COND^./V N-l M NOZ^Vn-1 


T = - 1 — 1 1 1 


N 


gN N n _ , . 

A k'V 1 ' 0 


c'N-1 ~ ^COND^N ^ M NOZ^N 


takes into account the change in gas temperature due to evaporation, 
condensation, and gas efflux from the nozzle. 

In this analysis, the gas temperature will be determined from: 


N n 


T = 
gN 


S 3VA-1 (M c>N-f ®0 OTD>N ( Vn- 1- (M NOZ>N ( VN-r < M s V T gVj 


^ G ij + ^ M c^N-l ^COND^N ^ M NOZ^N 

-e/r(t ).. . 

AHKC c m e g N 1 . h A (T -T) 

+ At + ^ — At 

pC N m C N 


(hA + WC ) (T-T ) 
P o 

m C 


h 0 A (T-T ) 

_ _ w . 2_w ... w 

N m C N 


where (M g )^ - mass of solid or liquid intermediate compound formed during 


th 


the N time interval. 


Pressurization as a function of time will then include the effects of 
evaporation, condensation, mass flow out the nozzle, intermediate compound 
formation, gaseous chemical reaction, and heat transfer to and from the wall. 
The pressurization equation then becomes: 


A Pi A 


ar Wl^ W J % 
^dhjN-l 
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a. A At 


/--M g 


c N 27 tR(T )„ . 

g N-l 


( Vn-1 A * 


y+1 


(YU ' 2 7 1 


‘7n 


g N-l 


V R y S V (M ! ) N lt R' T~ Lt^’ T 3 ": J 

gN gN-1 


Ignition process. - The transient model being developed will deter- 
mine the conditions in the chamber at all times up to ignition. 

The pressure generated by either a detonation or an explosion can be 
calculated. The equation P /P =(1+Q'y M 3 )/l +y can be applied to determine 

M 1 1 1 fa 

the magnitude of the transient overpressure due to a detonation. The use 
of this equation will depend upon whether the detonation induction distance 
and engine contraction ratio (A c /A t ) criteria are satisfied. 

The magnitude of a transient overpressure due to an explosion can be 
calculated from the ideal gas law. This can be done by assuming that the 
propellants react instantaneously throughout the entire chamber volume or 
that they collect and react in only a small segment of the chamber. 

Mathematically, this problem area has been reasonably solved. The 
magnitude of the detonation or explosion can be determined from conditions 
existing at ignition. The identification of the triggering mechanism will 
probably have to be determined experimentally. However, studies can be made 
of different spiking mechanisms (such as formation of hydrazinium nitrate, 
explosion or reaction of frozen propellants on the chamber wall, etc.), 
so as to help design an engine minimizing the occurrence of excessive 
initiating transient overpressure. 
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Results of Chamber Pressurization Computations 


By means of the chamber pressurization program, the gas pressure, 
gas temperature, mass of propellant evaporated, mass of gas condensed, 
mass efflux from the nozzle, liquid drop temperature, and fraction of a liquid 
drop frozen can be computed for the nitrogen tetroxide-Aerozlne 50 system. 

The applicable computer program listing is given in Appendix A. In its 
present state, the pressurization program treats each propellant separately, 
that is, pre-ignition chemical kinetics, the related chemical energy release 
and temperature change have not been incorporated in the computer code to 
date. However, these interaction effects have been formulated analytically. 
It should be emphasized that the chemical kinetic parameters for the impor- 
tant hypergolic propellant reactions are difficult to determine, but that 
meaningful parametric studies can be performed within the framework of the 
analytical model. 

Parametric studies. - Several parametric studies, the results of which 
are given graphically in Figures 13 through 16, were carried out utilizing the 
present pressurization program. In these studies, the dependence of the 
chamber pressure on the flow rate, chamber wall temperature, accommodation 
coefficient, and droplet radius is determined. 

Figure 13 shows the variations of chamber pressure history with 
flow rate for the oxidizer (^O^) and fuel (Aerozine-50) separately. A 
40 percent increase in flow rate from 0.345 lb/sec to 0.483 lb/sec brings 
about a maximum absolute increase in chamber pressure of «s0.33 psia for 
the oxidizer at«2.3 msec after injection. In the case of the fuel, the same 
change in flow rate yields an increase in chamber pressure of «0.13 psia 
which does not vary with time afters 2 msec. 

The calculations of the effect of wall temperature on chamber pressure 
(Figure 14) show that for N„0„ during the first msec, a relatively small 

L H 

difference in pressurization rate is obtained, whereas after «4 msec at least 
at 1 psia difference in chamber pressure is obtained for a 91°F change in wall 
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temperature. 

Figure 15 shows the significant effect of the accommodation coefficient 
on early pressurization rate and on the absolute pressure values obtained after 
a few msec following injection. This figure applies to N O and is based on 

c* “ 

an initial droplet temperature of 530°R. 

The effect of oxidizer droplet radius on chamber pressurization is 
demonstrated in Figure 16. As one would expect, smaller drop radii enhance 
chamber pressurization. Since the oxidizer flow rate is only 1/4 of the highest 
value employed in the calculations, both sets of droplet radii, even after 
7 msec, yield a comparatively low chamber pressure, that is, below 1.5 psla. 

The above computational studies have demonstrated the relative effect 
of several important parameters governing the physical processes taking place 
within the combustion chamber. Further development of the mathematical model 
will make possible the quantitative study of additional chemical and physical 
interactions, and thus help to define engine design and operating conditions 
minimizing or eliminating spiking. 
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CHAMBER PRESSURE PSIA CHAMBER PRESSURE PSIA 



MILLISECONDS 

FIGURE 13 . EFFECT OF FLOW RATE ON PRESSURIZATION. 



FIGURE 14 . EFFECT OF WALL TEMPERATURE ON OXIDIZER PRESSURIZATION. 
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CHAMBER PRESSURE PSIA _ CHAMBER PRESSURE PSIA 


ALPHA = 1.0 


R, = 83XI0"6fT. 

R 2 = 25X10~ 5 FT. 

R 3 = 72XI0“5 ft 



3.0 4.0 

MILLISECONDS 

FIGURE 15. EFFECT OF ACCOMMODATION COEFFICIENT ON OXIDIZER 
PRESSURIZATION . 



MILLISECONDS 

FIGURE 16. EFFECT OF DROPLET RADIUS ON OXIDIZER PRESSURIZATION. 
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VI. CONCLUSIONS 


The following enumeration presents significant conclusions which 
were drawn from the experimental, order of magnitude, literature survey, 
and computational work concerning overpressure phenomena carried out in 
this study: 

(1) Experiments showed that cavitation and two-phase flow is present 
within the injector during the space starting sequence of an engine, 
and that this phenomenon may significantly limit the initial flow 
rates. 

(2) It was shown that especially in regions close to the injector face 
high amplitude detonation waves can be supported by a two-phase 
propellant mixture, whereas one-phase, gaseous detonation 
produces relatively low detonation pressures anywhere in the 
chamber. Thus, in the nitrogen tetroxide-hydrazine system 
observed high amplitude overpressures may only be accounted 

for by two-phase detonations. 

(3) Heterogeneous detonation calculations show order of magnitude 
agreement with maximum spike pressures observed in certain 
engines . 

(4) Existing literature indicates that intermediate reaction products 
(gaseous, liquid, and solid) may play an important role in spiking. 

(5) It is not clear whether transient overpressures stem from material 
(intermediates or pure propellants) accumumulated on the walls or 
originate from reactions in the combustion chamber space. This 
points to the importance of understanding the triggering mechanism. 

(6) The mathematical model developed in this study, in its present 
state, is capable of predicting chamber conditions during the 
preignition time period. The computer solution may be used to 
perform parametric design studies. 
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APPENDIX A 


COMPUTER PROGRAM LISTING 


PROGRAM PRESS ( INPUT»QUlPi'T, TAPE h*INPUT,TAPE 6*0UTPUT> 

TRaNSIEN] PRESSURE HISTORY PROGRAM 

COMMON /A/ TFP (2 ) » T AU ( 2 ) , t AUF 12 ) »TAUS(2) ,XMC2) »GAM(2) ,TG(2> , 

;PG<2>» M( 2 ) ,GfcVAP(2» ,Xi (?) ,xR(4,?> . ROL (?) ,CPL (2) »TO (2) 

COMMON / .HCOtK/X (4»50n«2) , T ( A , 500 » 2 ) »P (4 .500 * 2 ) , R ( 4 , 5 0 0 * 2 ) 

0 1 MENS I ON PV(IO), TD(JO) *TAG(?) »PV- ( 2 ) * PVF (4 ) * TDK ( 4 ) 

DIMENSION M] (20*?) .T1 <?0»?> *W (500 *2) 

DIMENSION T A ( <♦ ) »CPLA(4) *M?> ,r,COND(2) ,XMN0Z(2> *TQ(2) 

DIMENSION PG l t 2 ) 

DIMENSION ROP(?) 

DATA (R0P(1)=^9.Q?) , (POP (?) =56, 3?) 

DATA ( (PV (I) » l *1,9) *389. 29, 720., 1324,8,21 28. 32 *2733.1 ?,*4 19. 36* 

1 556 1 .28,6946.56, 1 0673. ?8) 

DATA ((TO (I) ,I*l,9)a4n.«<nA9O.*SI0.,530.,540.,560.,570.,580,.6l0. 

I ' 

OAtA ( (PYF(I) ,1*1, 4) *79.?, 300. 96, 720.. 105 1.2) 

DA | A ( (TDK (I) .1*1 ,4)*474.,528.,564.,582.) 

DATA ( (CPLA(I) ,I*l,4)s.3*67,.3696,.3752*.3825) 

DATA ( (TA(I) »i *1,4 >*524. 8,530. ,540., 550.) 

1 RE 0(5,9000) NN,m(1) ,M(2) ,NST0P, NO- 

DO 2 L-'l.NN 
00 2 J* UNSTOP 
X ( L * U *0. 

’ ( L » J , 2 ) *0 * 

T (L,0*1)*0. 

T (L, J,2) *0. 

P<L,J»1>=0. 

P(L» J,2) *n, 

R (L, J,2)*0. 

2 «(L*J*l>aO, 

HEAD (5, 9u 05) n(1,1)»W(1«2) ,TFP (1 ) »TFP (2) , T AUF < D * 

1 TaUF (?) ,TaUSU) , T aUS ( 2 ) ,XM(1),XM(2) *6AMM ) ,GAM(2) , TG ( 1 ) , Tfi (2) . 

?PG ( 1 ) ,PG(2) »Pvrf<1 ) *PVw ( 2 ) »ASTAR,YC,AC,ALPHA,DFLTN, (Xo ( I , 1 ) , I*l,NN) 
3, (XR ( I ,2) ,1* I »NN) 

READ (5,9010) (R(I.1,1) ,1*1 ,NN) 

HFAD (5,9010) (R ( I « 1 ,2) « 1*1 ,NN) 

HEAL) (5.9010) (T (1.1,1) ,1*1 ,NN) 

HEAD (5,9010) (T (1*1,2) ,1*1 ,NN) 

HEAD (5,9005) (>,1(1,1) ,1 = 1 ,NOW) 

HEAD (5,9005) (#1(1,2) ,1*1, NO*) 

READ (6, 9 005) (T1 (1*1) ,1*1 »NOw) 

REALM 5, 9 005) ( T 1 ( I ,2) ,1*1 .NO- ) 

9000 FORMAT <6I6> 

9on*n F Ohm AT (ot j.2.6 ) 

9010 FORMAT (4E12. 6) 

c compote Initial conditions 

HC * 1545. 


A-l 



nnn ctnct 


GCONQ ( 1 ) *0 • 

GCONO { 2 ) *0 , 

CONST i * 4,*ALPHA«DELTN*3. 1415/SORT (6.283l853*RC/32.2> 
C0NST2* 3,*ALPHA*0ELTN/SQRT (6.2831853*RC/32.2> 

C0NST3* C0NST2*AC/3. 

C0NST4 * ASTAR*SQRT (32. 2/1545. )»0ELTN 
NPaGE«1 

timeuo. 

TIME*0. 

N22«N0#-1 

NSTOP * nST0P*1 
TAG(1)»0. 

Tag (2) *o . 

DT1000»OELTN*1000. 

SELECT THE CORRECT PROBLEM 

NPP a 0 
NPF « 0 

I F ( M ( 1 ) . EQ . 0 ) NPF * ? 

IF (M(2) .EQ.01 NPF * 1 
IF (NPF ,EO. 0) NPF*3 
IF (M(l) . EQ.M (2) ) NPF «4 

IF( (M(l) .LE.m( 2) ) ,ANO. (Md >.GE.1>> N«l 

IF t <M (2) ,LE.M(1) ) .AND, (M(2) .GE.l) ) N*2 

IF (NPF «LT .3 ) N»NPF 

THE TIME LOOP 

00 lUOO Lai. NSTOP 
Tag (N> *o • 

GEVAP(l) «o. 

‘ GEvAP (2) «0. 

Ki(l) «n. 

xl(2) *0 • 

L**L 

IF (NPF. NE. 4) GO TO 2P5 
50 N«1 

T4G(N)»0. 

L •L-M(i)*l 

I F < (LX. GE.l) .AND, (M ( 1 ) .GT.O) ) GO TO 205 
100 n«2 

TAG (N) fO , 

LX«L-M (2) *1 

IF( (LX.GE.l) ,?ND. (M<2) .GT.OM GO TO 205 
GO TO 800 
1*5 BBaSQRT (XM(N) ) 
no ?00*NNN«1 iLX 
J»L -NNN*1 
CCa* ( J*N) 

DO 700 lel.NN 
IF (N.EQ.l ) GO TO 250 
00 251 K«1 * 3 

IF((T(I»J.2).GE.T0K(K)).aND.(T(I.J»2).LE.TDK(K*1))) GO TO 252 
251 CONTINUE 
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T.TO 


2F»? P < \ ♦ J.2) * (PVF (K*l ) -PVF (k) )• (TOK < K ♦ 1 ) -T U • J*2) ) / < TOK ( K ) •* TDK t K. 1 ) ) ♦ 
1PVF t K * 1 ) 

TAU(2)s425.8 

!M1 (1.J.2) .LT.140.) ROlJZls 68 , 64*T ( I . J . 2 ) 

IF (T < I * J ♦ 2 ) . GE « 1 4 0 , ) ROL<?)* 03158*T ( I .J.2) *7.1.06 

CPL<2) * ,000?5*T (I. J.2)., 56 


c 


GO TO 31 

260 CONUMUt 

IFtTU.J* | ) ,LT, 471,84) GO TO 280 

UO 2bo K® ) . 8 

2*0 CONTINUE* * * #GfL,TD(K> T(I,j,i),LEJO(K*in GO To 270 

270 CPV(K».PV<R*1))«CTDfK*])-Tll f J,i))/ t TOIK*l)-rO(K)I*PVCK4 

GO 10 290 

2f)r) p U.J.l) S 6,l58*TU.J» 1J-2624.05 


2=>o 


TAU(1> * 178. 


TF( T (I » J* 1 > ,LT. 140,4) ROL(l)* 122,49 

TF ( (T ( I » J, I ) ,0E, 140.5) ,4Mn. (T (I * J» 1 ) .LT.471 .84) ) POL ( 1 ) *-,08433* 

IT <I»J.l)*l34,33 

IF (T (I#J. \) .GF. 47 1,84) ROL(1)s-,07934*T(I.J.1)*131,9p 
I F ( T ( I * J « 1 ) ,G£«524.<) GO TO 300 

( , Tt f ’ J * l) • GE ' 47 l* p ^ ) - ANO , ( T < 1 • J* l ) .LT.524.8) ) CPU1) a, 0002325* 

) T (T * J ♦ ). ) * , 2 4 4 b 

I* ! ( Vi ,J,1) • ge *1°5* £ 1) .and, (T ( If J.l) .LT.471,84) ) CPL (1 )«, 0004909 
1*T < I * J. 1 ) ♦, 0bb9 

IF '0 T 0 1 3 ^ 0 1 * • LT,105,91) CPLH)* .001?5*T (l.J.n-. 02349 
3o0 00 305 K r; 1 , 3 

305 !F( < C0NT1nUe' IJE ' TA<K)J * ANi °* .LE.TA(KM>> ) GO TO 307 

3o7 CPL ( 1 ) = (CPL A (K) »CPLA («;♦! ) )*(TA(K*1)*T(I.J.1))/(TA(K*1)«’TA(K))* 

1 CPLA(K*1) 

310 CONTINUE 

IF ( < ( I . J.N) ,GT, 1 , ) GO TO 312 

G(N) * UP ( 1 1 J.N) -PG (N) ) *R < I , J.m) **?/S0RT ( T ( I » J,N) ) ) « CONSTl 
0 ( N ) = G(N)«BH 
IF(G(N) ,LT.O.) G(N)sQ, 

GO TO 3 16 
312 G(M)xO, 

316 TOG * (G (N) *CC*XR( I *N) / (POP (N) *» ( I » l.N) #*3) ) #,p3«7 
GeuArvNi = GEvAPiNi* TOn 
TAG(N) s T4G (N) * T0G*T(T.JtN) 

IF ( X ( I » J.NJ..GT ,0 , ) GO TO 315 

IF (P (I * J.N) ,LE,PG (N) ) GU TO 315 
IF ( TFP (N) ,EQ. T ( I . J.N) ) GO TO 3)5 

«(T.J.1.\>) x W (T . J.N) * ( 1 08*{P(I»J.N)-PG(N) ) *rONST2/ (POL (N) 

1 *SQ«T ( f U.J.N) ) *9 ( X . J.N) ) >**,3333 
GO TO 3 2 1 ' 



»N) « R ( I « J»N) 

3?0 ; IF(.X(I*J»N> .GT.O.) GO TO 360 
' : fF(P(I*U*N) .GT.PG(N) ) GO TO 340 

, 1 # n > • t ( i # j ♦ n ) 

350 

3^e f<ItJ*liN)«T(I*J,N)-C0NST2*TAU(N)*(P(T.J.N).PG(N))/(R0L(N)*CPL(N)* 
V I * *SQRT ( T ( I ♦ J*N) > )*B8 

; vi; ; f?tTUfj4i*N>,i.T.TFP(Nn t<i*j*i*n)«tfp<n> 

X(!tJ«N) ♦ <ri(N)«TAUS<N)/(ROL (N)*R(I.J,N)**3*TAUF(N) ) ) 


B TO 370 

fJ#l»N)«TFp(N) 
JO 350 

CONTINUE 


CONTINUE 




K*1,N22 

(K,N> .AND. TIME l .LE. T1 <K*l *N> » GO To 510 


: ,T> • '• 


' ;#3| 


|pi|f#V«tWl(K*ltN)*Wl (K,N) )*(Tl (K*ltN) -TIMED /(Tl (K,N) -T1 <«*1#N 

|F (NPF.EQ. 1 ) GO TO 800 
IF(NPF.EQ.2> GO TO 800 
IF(NPF.EQ.3) GO TO 710 
JF( (NPF.EQ, 4) .AND. (N.FQ.Z)) GO TO 800 
, IF ( (NPF.EQ. ♦) .AND.(N,EQ.1D GO TO 100 
IF(N.EQ.l) MPF«2 
fe*r|F4Nri.Q«2) MPF*1 

|F (M(MPF ) .EQ ,U NPF ■ ' 
fep.iWF*NE,*) GO TO 800 
Wif .EQ.2) GO TO 100 

.TO 50 


(»*»• t (l*G (N> -PVW (N1 1 4C0NST3/SQRT (TG (N) D«BB 
(BCONO(N).LT.O.) GCOND <N)«0. 

- v^I^ fj#i01(NI m (PG (NUSQRT (TQ (N) D *C0N5T4#BB*SQRT (GAM (N, * (2./ (GAM (N) ♦ 

ilWPP # ^aF M < N >*lt)/<GAM(N)-l.)n 

*§'(*>'• ( TAG (N) ♦VC*XM(N)*PG(N)/RC -GCOND (N) «>TG IN > -XMNOZ (N) *TG ( 

(GEVAP (N) *VC*XM (N) OPG (N) / (RC *TG (N) ) -GCOND (N) .XMNOZ (N) ) 

* 1*GEVAP(N) «RC*TQ (N) / ( VCOXM (N) ) -RC*TQ (N) *GCOND (N) / ( VC*XM (N) ) 
fB(Nl*XMNOZ(N)/(VC*XM(N) ) *PG (N) *TQ (N1 /TG (N) 

(N)«PG(N)/144, 

'IF (NPF. NE, 4) GO TO 820 
• IF ( (NPF.EQ. 4> t AND. (N.EQ.2D GO TO 815 
N«2 

■ t 60 TO 810 

;fPGTOTL»PGl ( 1 ) ♦PGT (2) 

I GO TO 82^ 

#j|*##«1 , O.TL«PGl (N) 
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H?5 IF(NPAGE ,E0. 1) WRITE (6.2000) 

I F ( nPF ,FQ. 4) GO TO ono 

IF ( N ,EQ, 1) GO TO 920 

wR I re (6*2010) TIME »PG1 (2) »GEVAP(?) .GCOND (?) *TG(2) 

NPAGfcaNPaGE*! 

IF (NPAGE.GE.60 ) mPAGE*1 
GO TO 99 > 

9?o *RJTE (6,2020) T I m£ ,PG1 U ) .GEV *P <1 ) , GCOND « 1 ) * TG ( 1 ) 

NPAGE*NP«GE*1 
IF (NPAGE.GE.6-i) NPAGEsl 
GO TO 990 

9o0 wRl TE (6*2030) TIME.PG1 (1 ) ,PGl (2) » GEVAP ( 1 ) ,GEVAP<?) *QC0N0(1> tGCOND 
I (?) »TG(1) ,TG(?) 

MPAGEaNPAGE*! 

IF (NPAGE.GE.6-i) NPAG£ = 1 
9Q0 CONTINUE 

2000 F0RMAT(1H1.4A.4HTIME.9X.GHPGM ) ,8X,5HPG<2) ,8X,SHGE\/Ap (1 ) ,5X»8HGEVA 
iP<?) »5X.0HGCONO(1 ) » 5X • 8HGCONO ( 2 ) ,5X,5MTG(1) ,«X*5HTG<?> ) 

2010 FORMAT (IH .EU.4.I3X.F 13,4.1 3Xt El 3.4,1 3X.E 13.4. 13X,E 13, 4) 

2020 FORMAT ( 1H , E 1 3 .4 ,E 1 3 .4 , 13 X ,El 3 .4 , 1 3X ,F 1 3 .4 , 1 3X ,E 1 3 .4 ) 

2030 FORMAT 1 1 H ,9El3.4) 

3000 CONTINUE 

TG(l)sTQ(l) 


T G ( 2 ) * T Q ( 2 ) 
TIME«TIME*oT1000 
timei«timei*oeltn 
1000 CONTINUE 
GO TO i 
Eng 
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